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Implementation example: 

(read discretization time step) FIG.13A 
get dt 

[reset oscillators] 

130 for each N [loop over resonators] 

get omega(N) [read frequency *2*pi] 
get tau(N) [read time decay constant] 

get weight(N) [read weight] 

u 0 (N)= (2-(omega(N)dt) 2 ) / (1+dt/tau (N) ) 

[one example of equations describing 
[a damped harmonic resonator] 

u m (N)= -(1 - dt/tau (N) ) / (1+dt/tau (N) ) 

u.(N)= dt / (1+dt/tau (N))/tau(N) 

[u are constants] 

bm(N) =0 [b are amplitudes of resonators] 

bO(N) =0 [ bm is amplitude of previous time step] 

[bo is amplitude of actual time step] 

131 done with loop over N 
[loop over time steps] 



132 



for each t 

input ap [read input signal a of the] 

[next time step] 



if (first t) then [the first time, reset stack] 

[for input; otherwise a would] 
[not be defined] 

135 °o=ap 
am-ap 
end if 

[calculate difference aprime of input signal a] 
aprime - (ap-am) 

continue with FIG.13B 
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continuation of FIG.13A FIG.13B 
[propagate resonators; loop over all resonators] 
for each N 

bp = u 0 (N)*b 0 (N) + u m (N)*b m (N) + ui(N)*aprime 

[resonator for following time] 
[step] 

b m (N) = b 0 (N) [stack is updated/shifted] 

b 0 (N) = bp 
done with loop over N 

[at this point the weights weight (N) can be recalculated] 
[as a function of the b 0 (N) ] 

[In this implementation example the weights are constant] 

[weight and combine resonator output signals to generate] 
[output signal c] 

c = 0 
for each N 

c = c + weight(N) *b 0 (N) 
done with loop over N 

[reset stack for input] 
a m = a 0 
ao = ap 

[output output signal c] 

output c 
done with loop over t 

end. 



